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We theoretically find that in a multiferroic chiral magnet Cu 2 0 Se 03 , resonant magnetic excita¬ 
tions are coupled to collective oscillation of electric polarization, and thereby attain simultaneous 
activity to ac magnetic field and ac electric field. Because of interference between these magnetic and 
electric activation processes, this material hosts gigantic magnetochiral dichroism on microwaves, 
that is, the directional dichroism at gigahertz frequencies in Faraday geometry. The absorption 
intensity of microwave differs by as much as ~30% depending on whether its propagation direction 
is parallel or antiparallel to the external magnetic field. 

PACS numbers: 76.50.-|-g,78.20.Ls,78.20.Bh,78.70.Gq 


Collective excitations of spins in magnets, so-called 
magnons or spin waves, can be activated not only via 
a direct process with ac magnetic field H‘^ coupled to 
magnetizations but also via an electric excitation by ac 
electric field E‘^ coupled to charge degrees of freedom. 
When the magnon or spin-wave modes have simultane¬ 
ous activity to the and E‘^ components of electro¬ 
magnetic waves, interference between the two activation 
processes, that is, the magnetically activating and the 
electrically activating processes, gives rise to peculiar op¬ 
tical and/or microwave phenomena, so-called optical ME 
effect. One of the most important examples is the direc¬ 
tional dichroism, that is, oppositely propagating electro¬ 
magnetic waves exhibit different absorptions. 

Multiferroic materials with concurrent magnetic and 
ferroelectric orders EHi provide an opportunity to re¬ 
alize the electric-dipole active magnons (so-called elec- 
tromagnons) [SE3 , and thus the optical ME effect via 
the magnetoelectric coupling 14l4lfil|. Indeed observa¬ 
tions of the directional dichroism have been reported for 
several multiferroic materials such as Ba 2 CoGe 207 [13- 
[l9l| . i^MnOa (i? =rare-earth ions) [^, [^, and 
CuFei_a,Gaa;02 [Hj, in which nontrivial spin orders in¬ 
duce the ferroelectric polarization via the relativistic 
spin-orbit interaction. In these materials, the optical 
ME effect is observed at the electromagnon resonance 
frequencies in the terahertz (THz) regime. 

The directional dichroism is observed also at higher 
frequencies, i.e., x-ray and visible-light regimes in several 
polar magnets, which is caused by electron transitions 


among the spin-orbit multiplets |23l - t29l| . However, ob¬ 
servations of the effect at gigahertz (GHz) frequencies 
are quite limited and the effect observed so far is very 
tiny whose difference in absorption intensity is only 2.5% 
at most [s^, while the directional dichroism at GHz fre¬ 
quencies is anticipated for application to microwave de¬ 
vices [1^ . This is because most of the well-known multi¬ 
ferroic materials based on simple spiral or antiferromag¬ 
netic spin structures with short-period modulation tend 
to have relatively large spin-wave gaps of several meV, 
which inevitably results in rather high resonance frequen¬ 


cies in THz regime. To achieve the microwave ME effect, 
long-period magnetic textures with tiny spin-wave gap 
should be examined, and slowly modulating magnetic 
structures induced by the Dzyaloshinskii-Moriya inter¬ 
action (DMI) are promising for this purpose. Indeed, 
the nonreciprocal directional dichroism of microwaves in 
Voigt geometry was theoretically predicted for the DMI- 
induced skyrmion phase (sij. 

In this Letter, we theoretically predict that electro¬ 
magnon excitations in magnetically ordered phases of 
Cu2 0Se03 exhibits unprecedentedly large magnetochi¬ 
ral dichroism at GHz frequencies, that is, the microwave 
directional dichroism in Faraday geometry. In the pres¬ 
ence of H, the conical spin phase or the field-polarized 
ferriinagnetic phase emerges in the bulk samples depend¬ 
ing on the strength of H. When H is applied in a cer¬ 
tain direction, and a microwave is irradiated parallel or 
antiparallel to H, the absorption intensities for the oppo¬ 
sitely propagating microwaves differ by as much as 30%. 
Such a huge directional dichroism at microwave frequen¬ 
cies has never been observed in a single material. This 
effect is traced back to the resonantly enhanced magne¬ 
toelectric coupling, and therefore essentially distinct in 
microscopic mechanism from the traditional microwave 
non-reciprocal device based on microwave polarization, 
potentially leading to a unique microwave device. 

The crystal and magnetic structures of Gu2 0Se03 con¬ 
sist of a network of tetrahedra composed of four Cu^+ 
(5=1/2) ions at their apexes as shown in Figs. [T](a) 
and (b). In each tetrahedron, three-up and one-down 
collinear spin configuration is realized below Tc~58 
K [H, . This four-spin assembly can be regarded as 

a magnetic unit, and is described by a classical magneti¬ 
zation vector rrii whose norm m is unity. We employ a 
classical Heisenberg model on a cubic lattice to describe 
the magnetism in a bulk specimen of Cu2 0Se03 [35l - [33| . 
which contains the ferromagnetic-exchange interaction 
and the Dzyaloshinskii-Moriya interaction among the ef¬ 
fective magnetizations rrii and the Zeeman coupling to 
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FIG. 1; (color online), (a) Crystal structure of Cu 20 Se 03 . 
(b) Magnetic structure of Cu 2 0 Se 03 . (c) Phase diagram 

of the spin model given by Eq. m with J=1 meV and 
_D/J=0.09. (d) Schematic figure of the conical spin struc¬ 

ture. (e) Symmetry axes in the chiral cubic P2i3 crystal, (f) 
Ferroelectric polarization P || [001] under H jj [110]. (g) Ab¬ 
sence of P under H || [010]. (h) Calculated net polarization 
Pc in the field-polarized ferromagnetic state under HPc as a 
function of the angle 6 between H and the h axis. 


the external H. The Hamiltonian is given by, 

7^0 = —J 'TTii • fnj — D rrii x • 7 

<i,j> 

- ( 1 ) 

i 

where 5 = 2 , and 7 runs over d, b, and c in the cu¬ 
bic setting. Here details of the exchange interactions 
within each tetrahedron are neglected and only the ef¬ 
fective interactions among these magnetic units rrii are 
taken into account. Such a coarse graining is justified as 
far as low-energy excitations with long-wave-length mag¬ 
netic modulations are considered. We take J=1 meV 
and £>/J=0.09 so as to reproduce the experimentally ob¬ 
served Tc and periodicity (~ 50 nm) of the conical state. 

FigureHJc) shows a phase diagram of this spin model at 
T =0 as a function of the magnetic field H, which exhibits 
a phase transition between the conical phase and the fer¬ 
romagnetic phase at 77=691 Oe, in agreement with the 
experiments for the bulk specimen [381 [39j. Schematic 
figure of the spin structure in the conical phase is shown 
in Fig.[TJd), in which both the propagation vector Q and 


the net magnetization M are parallel to H. Note that 
the spin textures considered here are slowly varying in 
space, and thus the coupling to the background crystal 
structure is significantly weak. This justifies our treat¬ 
ment with a spin model on the cubic lattice after the 
coarse graining of magnetizations. 

The presence or absence of the ferroelectric polariza¬ 
tion P and, if any, its direction can be known from the 
symmetry consideration 38, 4^. The crystal structure of 
Cu 20 Se 03 belongs to the chiral cubic P2i3 point group, 
which has three-fold rotation axes, 3, along (111), and 
two-fold screw axis, 2 i, along ( 100 ) as shown in Fig.[T][e). 
This crystal symmetry is not polar, and thus there exists 
no spontaneous P. Although the conical and the ferro¬ 
magnetic spin states are not polar, either, combination 
of the crystal and the magnetic symmetries renders the 
system polar, and allows the emergence of P. As shown 
in Fig. [T][f), the emergence of P||[001] perpendicular to 
the net magnetization M{\\H) is expected for the coni¬ 
cal and the ferromagnetic states formed under ff||[ 110 ] 
since only the 2 i axis along [ 001 ] remains as a symmetry 
axis. On the other hand, the emergence of P is forbid¬ 
den under ff||[ 010 ] since three 2 i axes survive as shown 
in Fig.HDg). 

Microscopically the local polarization pi at the ith 
tetrahedron is given using the magnetization components 
rriia, rriih, and as. 


Pi — i^Pia: Pib^ Pic) — ^{pflih^^ici'bblic‘bbliaT'bblia^^ib) • ( 2 ) 

The net magnetization M and the ferroelectric polariza¬ 
tion P are calculated by sums of the local contributions 
as and respectively, 

where the index i runs over the Cu-ion tetrahedra, N is 
the number of the tetrahedra, and H(=1.76x 10“^® m®) is 
the volume per tetrahedron. The coupling constant A is 
evaluated as A=5.64 x 10“^^ /iCm from the experimental 
data [3l|. Figure |T][h) shows calculated net polarization 
P in the field-polarized ferromagnetic state when H is 
applied within the c-plane as a function of the angle 6 
between H and the b axis (see the inset). We find that 
the positive (negative) P emerges along [001] when 9 is 
negative (positive). 

According to Fig.[T][h), one realizes that oscillation of 
M (AAi’‘^||[100]) induces oscillation of P (AP“||[001]) 
when AT||iT||[010] [see Figs. Ill[a)-(c)], and conversely 
AP“||[001] induces AAf^HjlOO] via the magnetoelectric 
coupling. This means that both fi’‘^||[100] and P‘^||[001] 
components of microwave can activate the coupled oscil¬ 
lation of M and P (Note that the response time of M 
(P) against the change of P (M) is governed by the elec¬ 
tron transitions among the orbital multiplets, and thus 
is much shorter than the typical time scale of the oscil¬ 
lations). To see this, we calculate the following dynam¬ 
ical susceptibilities by numerically solving the Landau- 
Lifshitz-Gilbert (LLG) equation using the fourth-order 






































3 





_I_I_I_I_I_ 

0 1 2 3 4 5 6 

Frequency (GHz) 


FIG. 3: (color online). Imaginary parts of the calculated 
dynamical magnetic, dielectric, and magnetoelectric suscep¬ 
tibilities, Imx™, Imy®®, and Imy;J™, as functions of the fre¬ 
quency for the conical state at Jf=518.4 Oe and the ferro¬ 
magnetic state at 7F=1382.4 Oe when the static magnetic 
field ilf||[010](||b) is applied. 


and 


FIG. 2: (color online), (a)-(c) In the presence of net mag¬ 
netization M\\H under ilf||[010], oscillating magnetization 
component AAF^dKlOO]) is accompanied by the oscillating 
polarization component AP“(||[001]). (d) Conhguration of 
microwave and components, for which the magne- 
tochiral dichroism occurs under JF||[010]: fc“||±Ff, Ff“||[100] 
and F;“||[001]. 


Runge-Kutta method: 

AM‘^ 

magnetic susceptibility, 
dielectric susceptibility, 

AP^ 

= A ~ magnetoelectric susceptibility. 




eo 


electromagnetic susceptibility. 


Because of the symmetry, the relation X®'^(w) = 
holds. The LLG equation is given by 


dtrii 

dt 


—m, X HI 


eff 


OG 


m, X 


driXi 

dt ’ 




(3) 


where aG(=0.04) is the Gilbert-damping coefficient. The 
effective field JTf® is calculated from the Hamiltonian 
'H—'Ho+'H'it) as = —d'H/drrii. The first term "Hq 
is the model Hamiltonian given by Eq. ©• The pertur¬ 
bation term 'H'(t) represents a short rectangular pulse of 
magnetic field AH{t) or electric field AE{t), which are 
given, respectively, by. 


TL'{t) =-gfj.BtJ.o'P^AH{t) ■ (4) 

i 


n'{t) = -Y,AE{t)-p,. 


(5) 


After applying the pulse at t=0, we trace time profiles 
of M{t) and P(t), and obtain their Fourier transforms 
AM„ and AP“. Dividing these quantities by Fourier 
components of the field pulses, AH‘^ and AE‘^, we ob¬ 
tain the susceptibilities. The calculations are performed 
using a system of 7V=20x 20x140 sites with the periodic 
boundary condition. 

In Fig[3l we display imaginary parts of the calculated 
dynamical magnetic, dielectric, and magnetoelectric sus¬ 
ceptibilities, i.e., Imy™, Imy®®, and Imy®™, as functions 
of the frequency for the conical state and the ferromag¬ 
netic state when the static magnetic field ff|| [010](||b) 
is applied. These susceptibilities have resonant peaks in 
the GHz re gim e in agreement with the microwave exper¬ 
iments 411. Moreover all of these dynamical suscep¬ 

tibilities in each phase have peaks at the same frequency 
indicating that the resonant modes have simultaneous 
activity to H'“||[100](||a) and F;“||[001](||c). 

This magnetoelectric activity causes the microwave 
magnetochiral effect in Cu 20 Se 03 [see Fig. [^Kd)]. For the 
electromagnetic wave, the relation H^\\k‘^ xE‘^ holds, 
indicating that relative directions of the H‘^ and E‘^ 
components are determined by the propagation vector 

, and their relationship should be reversed upon the 
sign reversal of k‘^. When the lineally polarized elec¬ 
tromagnetic wave with H^\\a and £“^110 propagates 
parallel (antiparallel) to H\\h, that is, sgn[Refc‘^]=-|-l 
(sgn[Refc‘^]=—1) with fc“=fc“b, the and £“ com¬ 
ponents contribute in an additive (a subtractive) way 
to excite the coupled oscillation of M and P, resulting 
in weaker (stronger) absorption of the electromagnetic 
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FIG. 4: (color online). Calculated frequency-dependence of 
the absorption coefficients, q+ and a_, for microwaves with 
sgn(Refc“)=-Fl and sgn(Refc‘^) = —1, respectively, at several 
values of H when A:‘^||iF||[010], iF‘^||[100] and f;“||[001] in 
the conical state (a) and the ferromagnetic state (b). 


Magnetic Field (A/m) 



FIG. 5: (color online). Galculated magnetic-field dependence 
of resonance frequency and magnitude of the magnetochiral 
dichroism Aa/o-. 


wave. 


The expression of the complex refractive index N{ui) 
is derive d by solving the Fourier-formed Maxwell’s equa¬ 
tions as [17|, 


N{uj) = -k‘^ 

UJ 


- V'[ecc°“+X“(w)][/iaa“ 

-F sgn{Rek‘^)[Xac{^)+X 


xsrM] 
»]/ 2 , ( 6 ) 


for fc‘^=fc“6||[0I0], iT‘^||[100]||a and £;‘^||[00I]||c. The 


absorption coefficient a{uj) is related to N{uj) as. 


, , 2ujk(u!) 
a{uj) = - 


2uj 


ImA^(a;), 


(7) 


and thus attains the directional dependence via the sign 
of Refc“. Here k{lu) = ImiV(w) is the extinction coeffi¬ 
cient. 

Figures lUJa) and (b) display calculated frequency- 
dependence of a+ and a_ in the conical phase and the 
ferromagnetic phase, respectively, for several values of 
H, where and a_ are the absorption coefficients 
for microwaves propagating parallel and antiparallel to 
Ff||[010], respectively. In the calculation, we take e^=8 


and according to the experimental data 14, 

We find that the microwave absorption is resonantly en¬ 
hanced at the eigen-frequency of the electromagnon exci¬ 
tation, and there exists significant difference between 
and a_. 


Calculated magnitude of the directional dichroism 
Aa/a_ where Aa = a+ — a_ is plotted in Fig. [SJ This 
quantity is governed by the amplitude of ferroelectric po¬ 
larization P, which is nearly proportional to the square 
of the net magnetization M^. The dichroism increases 
in the conical state with the growth of P and M as the 
magnetic field increases. In contrast, once the system 
enters the ferromagnetic state, the saturated M gives 
nearly constant directional dichroism although there still 
exists slight field-dependence due to the field-dependent 
resonant frequency via uj in Eq. ©• The directional 
dichroism is enhanced at the phase boundary between 
the conical and the ferromagnetic phases, and its mag¬ 
nitude reaches as much as ~30%. Calculated electro¬ 
magnon resonance frequency plotted in Fig. O shows de¬ 
creasing behavior in the conical state, while increasing 
behavior in the ferromagnetic state as the magnetic field 
increases. This behavior is in good agreement with the 
experimental observations [321141| and coincides with an¬ 
alytical formula of the spin-wave gap given in Ref. g. 

In summary, we have theoretically predicted that 
magnetically ordered phases in the chiral multifer- 
roics Cu 20 Se 03 host gigantic microwave magnetochi¬ 
ral dichroism. This phenomenon results from interfer¬ 
ence between the magnetic and electric activation pro¬ 
cesses of electromagnons with GHz resonance frequen¬ 
cies. It has been demonstrated that long-period mag¬ 
netic structures in the chiral multiferroics without inver¬ 
sion symmetry can host gigantic dynamical magnetoelec¬ 
tric phenomena at GHz regime. In order to further en¬ 
hance the effect, search for novel chiral multiferroics with 
larger P or stronger magnetoelectric coupling is needed. 
For example, chiral multiferroics based on the inverse 
Dzyaloshinskii-Moriya mechanism Q as an origin of its 
P is worth trying to search because this mechanism tends 
to induce large P relative to the spin-dependent metal- 
ligand hybridization mechanism in Gu 2 0 Se 03 . Metama¬ 
terials, thin films, and synthetic nanomaterials are an- 
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other promising candidate to realize enhanced microwave 
ME effects through artificially designing large magneto¬ 
electric susceptibilities and intense electromagnon reso¬ 
nances [iH - lsoi . 
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